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Abstract 



We present theoretical results for (7, production on nucleons 

^ ' and nuclei in the kinematical region where the scalar isoscalar vrvr am- 

. plitude is influenced by the a pole. The final state interaction of the 

pions modified by the nuclear medium produces a spectacular shift of 
strength of the two pion invariant mass distribution induced by the 
. moving of the a pole to lower masses and widths as the nuclear density 

^ \ increases. 

The nature of the a meson has been the subject of continuous debate. Its 
^ , nature as an ordinary qq meson or a vrvr resonance has centered most of the 

^ \ discussion The advent of xPT has brought new light into this problem 

C ■ and soon it was suggested |^ that the a could not qualify as a genuine 

meson which would survive in the limit of large A*";,. The reason is that the 
^ , vrvr interaction in s-wave in the isoscalar sector is strong enough to generate a 

^ \ resonance through multiple scattering of the pions. This seems to be the case, 

and even in models starting with a seed of qq states, the incorporation of the tttt 
channels in a unitary approach leads to a large dressing by a pion cloud which 
makes negligible the effects of the original qq seed |Q . This idea has been made 
more quantitative through the introduction of the unitary extensions of xPT 
{UxPT) [^, 1^, 1^, P). These works implement unitarity in coupled channels in an 
exact form and use the input of the lowest and second order chiral Lagrangians 
of . The inverse amplitude method is used in [0 and an expansion in powers 
of O(p^) is done for the real part of the inverse of the scattering amplitude, 
while in the dynamics of the second order chiral Lagrangian is introduced 
via the explicit use of the exchange of genuine mesons, following the lines of 
Tm , and unitarizing with the N/D method. These works also justified the 



success in the scalar sector of the Bethe Salpeter approach used in 0, |T^. In 
all these cases the a meson appears as a pole of the vrvr scattering amplitude 
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in the second Riemann sheet, even when the second order Lagrangian, which 



contains information of the exchange of genuine mesons according to ||10|| , 
is neglected. These unitary models have been tested with success in many 
elementary reactions [p!3| . 

Another point of interest which can help us to understand the nature of 
the a meson is the modification of its properties at finite nuclear density. The 
importance of the medium modification of the vrvr interaction in the scalar 
sector was suggested in where the vrvr amplitude in the medium developed 
large peaks below the two pion threshold, somehow indicating that the a pole 
had moved to much lower energies. The issue has been revised and the models 



have been polished incorporating chiral constraints |T5|, [T^, with the result 
that the peaks disappear at normal density, but still much strength is shifted 
to low energies. 

Since present theoretical calculations agree on a sizeable modification in 
the nuclear medium of the tttt scattering in the a region, our purpose here 
is to find out its possible experimental signature in a very suited process like 
the (7, 7r°7r°) reaction in nuclei. Recent experiments at Mainz []TB[, where 



preliminary results show a very clear shift of strength of the invariant mass 
distribution of the two pions towards low invariant masses, seem to indicate 
that medium effects are indeed large. 

In the present paper we shall show how the tttt interaction in the final 
state of the (7, tt^tt^) reaction on nucleons enhances the cross section of this 
reaction and how the medium corrections on the tttt interaction in nuclei lead 
to an appreciable shift of the strength of the invariant mass distribution to- 
wards lower invariant masses. Although a similar shift has been claimed in 
the (vr, vrvr) reaction in nuclei |T^, |2D|, |2T[, the fact remains that there are still 
some discrepancies between these experiments and that the theoretical cal- 
culations [^, 23| do not reproduce the data. The reason is that the (vr, vrvr) 
reaction, involving initial pions which are much distorted in the nucleus, is 
quite peripheral and the effective densities tested are small. A possible way 
out to reconcile theory and experiment was suggested in |^ , showing that the 
small cross section in vr^p vr+vr^n at small invariant masses was abnormally 
small because of a subtle cancellation of large terms. A medium modification 
of these terms, through changes in the A^* properties and others not having 
to do with the vrvr interaction, could offset that cancellation and lead to larger 
final results. The (7, vr'^vr^) reaction is much better suited to investigate the 
modification of the vrvr interaction in the medium because the photons are not 
distorted and one can test larger nuclear densities. 

For the model of the elementary (7, vrvr) reaction we follow p4| which con- 
siders the coupling of the photons to mesons, nucleons, and the resonances 
A(1232), A^*(1440), A^*(1520) and A(1700). In the region of relevance to the 
present work, = 400 — 460 MeV, apart from some background terms, the 
A Kroll Ruderman term, diagram i) of Fig. 1 in p4[, is of importance and will 
be dealt with separately from the rest. The model of relies upon tree level 
diagrams. Final state interaction of the vrA^ system is accounted for by means 
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of the explicit use of resonances with their widths. However, since we do not 
include explicitly the a resonance, the final state interaction of the two pions 
has to be implemented to generate it. 

The 7iV — >■ Ntt^tt^ amplitude can be decomposed in two parts: the one 
that has as final state the combination of pions with isospin 1=0, first term of 
the RHS of Eq. ([I|), and the 1=2 combination, last two terms of Eq. 



|7r°(l)7r0(2) >= i|7r°(l)7r°(2) + 7r+(l)7r-(2) + 7r-(l)7r+(2) > 
3 

1=0 part 

^|7r°(l)7r°(2) + 7r+(l)7r-(2) + 7r-(l)7r+(2) > +|7r°(l)7r°(2) > 



1=2 part 

The interaction of pions in 1=2 in s-wave at these energies is very weak and 
hence we do not modify this part of the jN —>■ Nn^n^ amplitude due to the 
final state interaction of the pions. However, the 1=0 part is strongly modified. 
We have also checked that at the low energies involved here the pions come 
essentially in s-wave. 

In ref. [^] the renormalization of the / = (vr, vrvr) amplitude was done 
by factorizing the on shell tree level vrA^ — > vrvrA^ and tttt —>■ tttt amplitudes in 
the loop functions. This was justified in |]ll| for the tttt interaction. The same 
approach would lead in our case to 



1 (^^^0^0) 



(2) 



where Gt^tt is the loop function of the two pion propagators, which appears in 
the Bethe Salpeter equation, and t^^'' is the tttt scattering matrix in isospin 
1=0. In order to show clearly how one is lead to this equation, we show in 
Fig. the diagrams involved in the two pion production including their final 
state interaction. 






Figure 1: Diagrammatic series for pion final state interaction in 1=0 



The two final pions undergo multiple scattering which can be accounted 
for by means of the Bethe Salpeter equation, 



t = V + VG^J (3) 

where, following [|ll|], V is given by the lowest order chiral amplitude for 
TTTT — i> TTTT in / = and Ctttt, the loop function of the two pion propagators. 



3 



is regularized by means of a cut off in |11], or alternatively with dimensional 
regularization in [^]. In both approaches it was shown that V factorizes with 
its on shell value in the Bethe-Salpeter equation of Fig. |I|. Hence, in the Bethe- 
Salpeter equation the integral involving Vt and the product of the two pion 
propagators affects only these latter two, since V and t factorize outside the 
integral, thus leading to Eq. where VGT,T,t is the algebraic product of V, the 
loop function of the two propagators, Gtttt, and the t matrix. Coming back to 
Fig. |T] it is now clear that in case the vertex from where the two pions emerge 
is a contact term with a constant amplitude, the series implicit in Fig. ^ is 
summed as 



T + TG^^V + TG^^VG^^V + . . 
= T[1 + GUV + VG^^V + ...)] 



T(l + G^^t 

71 



(4) 



which is Eq. (|^). Now in the model for (7,27r) of there are indeed contact 
terms as implied before, as well as other terms involving intermediate nucleon 
states or resonances. In this latter case the first loop function in the diagrams 
of Fig. |1| is more complicated since it involves three propagators. Yet, if the 
intermediate baryon is far off shell, as is the case for most diagrams, then 
the baryon propagator does not change much in the loop function and the 
factorization of Eq. still holds. There is, however, an exception in the A 
Kroll Ruderman term, since as we increase the photon energy we get closer 
to the A pole. For this reason, and because its weight is important at these 
energies in the (7, vr+vr") amplitude which is needed in Eq. (^, we have singled 
out this term which we work out in detail below. The 1=0 part of the amplitude 
requires the combination of the diagrams shown in Fig. ^ a) b) 



A(1232) 



Jl+ P5 



Jl" P4 



+ A(1232) 



, Ji+ Pa 



TZ- P5 




a) 




b) 




71° P5 



71° P4 



c) 



Figure 2: Diagrammatic series for the A Kroll Ruderman term with final pions in 
1=0. 

and hence the vertex contribution is (assuming the A propagator, G^, the 
same in both cases) 



Ga{Sp5 S^e + Sp4 S^e) = GaS{pa + P5) S^e = GaSP S^e 



(5) 



where is the spin transition operator from spin 1/2 to 3/2 and e is the 
photon polarization. 

In the loop function originated from the sum of the terms a) and b) de- 
picted in Fig. 1^, assuming also the A propagator constant, we would have the 
contribution 
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^/ (277)4 2 5( 2 + g) Sh GA(lp^^p^_^2+j, (lp_g)2_m2+i, i^TTTT - 

= Ga 5P5^te'G^,(s) t^^(s) (6) 

since the term proportional to Sq vanishes because of parity reasons, and 
then we see exphcitly the factorization of the tree level amplitude. If we keep 
explicitly the A propagator in the loop some corrections arise since now we 
have the loop function with the A propagator and two pion propagators. This 
loop integral is performed following the steps of where the same loop 
with three propagators, albeit with only strong vertices, is evaluated. By 
performing explicitly the g° integration in Eq. (P) we obtain 



d\ ^ ^+^Ma 1 1 



(27r)3 '2 Ea 2u;u' + Lu + Lu' - LU - iu' + le 

1 1 



uJ + uJ'){^/s-EA-uJ-uJ')-uJP^] (7) 



where Ea = vMp + ^, lj = ^ml + (f + qf, uj' = ^ml + - g )^ 
Pa = f — g, Pa = p°+<^+'^ and ^/s is the CM energy of the initial photon and 
nucleon. 

The improved loop calculation results in a 10 per cent reduction in the cross 
section. The pion pole term (diagram j of Fig. 1 of ref. is essential in 
the (7,vr7r) model to guarantee gauge invariance, but is numerically negligible 
at the energies which we have here since it is proportional to P^, while the 
A Kroll-Ruderman term i) is proportional to P. Yet, inside the loops, the 
contribution can be bigger since the jmr vertex now is proportional to the loop 
momentum and not to the external pion momentum. In addition there is an 
extra p-wave ttA^A vertex which also goes like the loop momentum. However, 
there is also an extra pion propagator and one also finds cancellations from 
the poles of the different propagators. This was done explicitly in |^ where 



moderate effects from this pion pole term were obtained. We have performed 
the appropriate calculation and found that the loop diagram from the pion pole 
term is of the order of 10 per cent of the corresponding loop with the A Kroll- 
Ruderman term. Hence, we simply take care of it through the factorization 
approximation. We also take advantage of this numerical calculation to include 
two extra baryon form factors in the loop, as done in |2^, to account for the 



TriVA vertex correction. This is also done for the other diagrams since ttNN 
vertices are also involved. The inclusion of the form factors leads to a further 
reduction of about 15 per cent in the final cross section. Hence, the actual 
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Figure 3: Selfenergy diagram for the evaluation of the cross section including the 
final pious rescattering. 



modification of the amplitude should be 



T+(T' 



00 rpOO \ J. I ^ J. 



where T is the full (7, 7r°7r°) amplitude, is the (/ = 0, Is = 0) amplitude 
for (7,7r°7r'^), T'^j^ is the same as but calculated only with the A Kroll 
Ruderman term, Gtttt is the two pions loop including the strong irBB' form 
factors, t.jrTT is the vrvr scattering matrix in isospin 1=0, Gkr is the loop of Fig. |^ 
including in the integrand the two form factors and the T^j^ that depends on 
the momentum in the loops. (This is, Eq. (|^) including the form factors and 
the coupling constants and isospin coefficients of p4| ). 

After all this is done, we find it technically useful, in order to account for 
these elaborate loop corrections, to still apply the factorization of the (7, 27r) 
tree level amplitude of Eq. (0) but with a slightly modified form factor included 
in the (7,^77 loop function. The cut off of the monopole form factor is changed 
from 1 GeV to 625 MeV to implement these changes, including also the small 
effects from the projection in the two pion s-wave channel. This procedure is 
quite accurate numerically and prevents the numerical task from blowing up 
when we perform the calculations in nuclei. 

There is also a small technical detail. One of the terms in our approach 
contains the Roper excitation and its posterior decay into two pions in s-wave. 



The Lagrangian used is given in and the coupling constant is renormalized, 
as in |]23|, in such a way that when the vrvr final state interaction is taken into 
account the empirical N* width is obtained. 

The cross section for the nuclear process can be calculated using many 



body techniques in a similar way to |^, The method proceeds in two 

steps. In the first one the probability per unit length of a (7,277) process in 
nuclear matter is evaluated from the imaginary part of the photon selfenergy 
diagram of Fig. ^ corresponding to the cut of the horizontal line. By using 
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the local density approximation and, hence, substituting the density of nu- 
clear matter by the empirical density of the nucleus p(r) at a point r, this 
provides the reaction probability in this point of the nucleus. The second step 
requires to follow the individual particles produced through the nucleus in an 
inclusive process where the nucleus can be excited to any state. This second 
step is done by using semiclassical methods in which the pions produced follow 
classical trajectories and are allowed to undergo quasielastic collisions or pion 
absorption according to probabilities calculated previously using many body 
techniques. As shown in ||2^, this procedure was very successful in describing 
pion nucleus phenomenology at intermediate energies. 

With these ingredients the nuclear cross section is given by 



TT /■ 3^ /■ d^p r d^qi r d^q2 1 1 
kJ J (2^ J (2^ J (27r)3 2cj(gl) 2cj(gl) 

Si,Sf pol 

■ + E{p} - u;{q[) - ^(gl) - E{k + p- ql - (B)) 

■ F,{f,q-l) F2(r,g^2) (9) 

where n{p) is the occupation number for a density p(r). The factors Fi{r,qi) 
take into account the distortion of the final pions in their way out through the 
nucleus and are given by 



Fiif.qi 



exp 



roo I 

/ dli — /mn(rj 

Jr Qi 



(10) 



r + k qi/ \qi 



where 11 is the pion selfenergy, taken from where the interaction of low 
energy pions with nuclear matter was studied. This potential has been tested 
against the different pionic reaction cross sections, elastic, quasielastic and ab- 
sorption. The imaginary part of the potential is split into a part that accounts 
for the probability of quasielastic collisions and another one which accounts 
for the pion absorption. As we shall see, the probability that there is loss 
of pion flux through pion absorption at low energies is larger than through 
quasielastic collisions. One of the reasons is the Pauh blocking of the occupied 
states. 

When we renormalize the 1=0 amplitude to account for the pion final state 
interaction, we change G^^r and by their corresponding results in nuclear 
matter evaluated at the local density of the point r in the integral of 
Eq. (^. We do not include here the direct coupling of the two pions to the 
N* {14:4Q)h which was found extremely small in |]5D|. We do not include either 
the direct coupling of the two pions to a p — h excitation. The s-wave would 
proceed only through the tiny isoscalar vrA^ — ^ ttN interaction, and the p-wave 
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Figure 4: Contributions of the different isospin channels to the 27r mass distribution. 



part, weak in principle at the relatively small momenta of the pions, would be 
only relevant at excitation energies of about 30-50 MeV, far below the region 
of interest to us in the present problem. On the other hand, one should note 
that other alternatives have been proposed to study the modification of the 
isoscalar vrvr interaction in the nuclear medium based on the reduction of the 
pion decay coupling constant in the medium, which is tied to the dropping 
of the quark condensate via the Gell-Mann-Oakes-Renner relation pT| , 
or a combination of this effect and the dressing of the pion in the medium 
However, in the approach which we follow, one must be cautious not to 



include the dropping of / on top of the many body corrections done. Indeed, 
the change of / as a function of density in |4| and related works, obtained 
through the GOR relation, is linked to the renormalization of the time com- 
ponent of the axial current in nuclei, but in standard many body theory these 
currents are renormalized using the same lagrangians as in the evaluation of 



the pion selfenergy ^ Therefore, one should not modify / in 

these latter approaches to avoid double counting if using an explicit perturba- 
tive calculation with effective chiral lagrangians. Actually, a large change of 
/ at normal nuclear density would be difficult to accommodate to the quite 
well known pion nucleus phenomenology (pionic atoms, pion absorption, etc). 
Concerning this point it is worth noting the discussion in |^9|, ^ about the 
difference between the /(p) constant related to the quark condensate and the 
/ used in perturbative calculations with chiral lagrangians. These differences 
were also stressed in where it was also shown that taking the /(p) cou- 
pling related to the quark condensate in the many body evaluation of the pion 
nucleus optical potential led to unacceptably large widths of pionic atoms. 

We have also used the complex A selfenergy from |^ to dress the A 
propagator. In addition to the proper real part of the selfenergy in 



we 



add the effective contribution to the selfenergy 4:/9{f* / n^g'p coming from the 
iterated Ah excitation driven by de Landau Migdal interaction [HI, H3 . 
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[MeV] 



460 



Figure 5: Total cross section for tt^tt^p witli and without pion final state 

interaction. Experimental data from ref. 



In Fig. I we show the results of the invariant mass distribution of the 
two pions for the 7p — > n^n^p reaction. We can see the contribution of the 
1=0 part alone (TOO), the 1=2 part alone (T20) and the coherent sum of the 
two, both with and without renormalization (FSI) of the 1=0 amplitude. The 
renormalization of the 1=0 amplitude has important effects, nearly doubling 
the cross section. This is reminiscent of the similar enhancement found from 
chiral loops at threshold in When we sum coherently the 1=0 and 1=2 
amplitudes the shape of the distribution exhibits a double hump, one at low 
invariant masses and the other one at the high mass part of the spectrum. 
This shape is corroborated by the preliminary experimental results of [0 . 



The integrated cross section compared with experiment can be seen in 
Fig. 1^ where we can appreciate that the inclusion of final state interaction 
improves the agreement with the data. The amplitude for the (7, 7r°7r°) reac- 
tion on neutrons is calculated along the same lines and leads to a similar mass 
distribution as that of Fig. ^ albeit with a little smaller cross section. 

In Fig. ^ a) and b) we show the cross section for and ^osp^ assuming 
the FSI of the pions in the free case (p = 0), and we show the results without 
pion absorption or quasielastic steps, with pion absorption alone, and a third 
curve which corresponds to the case where the pions which undergo quasielastic 
collisions together with those absorbed are eliminated from the outgoing pion 
flux. This is actually not the case but the comparison of the two lower curves 
gives us a measure of the amount of quasielastic collisions undergone by the 
pions in their way out from the production point. The figures show that 
more pions undergo absorption than quasielastic collisions. The larger part of 
the quasielastic collisions would not change the charge of the pions, only their 
energy and momentum would be changed. In this case the two 7r° would still be 
there and their invariant mass would be somewhat changed, sometimes leading 
to larger two pion invariant masses and other to smaller ones. Hence, as an 
average the invariant mass distribution of the two pions should not be much 
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Figure 6: Two pion invariant mass distribution for 2tt^ photoproduction in ^^C and 
'^^^Ph. All three curves are calculated using the 27r final state interaction at density 
but they differ in the final pion distortion: continuous line: without absorption nor 
quasielastic scattering. Long dashed line: only final pions absorption. Short dashed 
line: final pions absorption and quasielastic scattering. 



modified by these quasielastic collisions. In other cases there could be change 
of charge and then we would not have two 7r° in the final state. However, this 
could also be compensated by having originally a vr+vr^ production followed by 
a collision of the 7r+ with charge exchange. The distortion of the final pions 
has been done simply by using the imaginary part of the optical potential and 
leads to a reduction of the cross section. The results for the case of removal of 
only the pions which are absorbed are obtained by putting in the imaginary 
part of the pion selfenergy, 11, in Eq. (0), the part which comes from the 
absorption and omitting the one that comes from quasielastic, which have 
been separated in Additional effects from the real part of the potential 
can be taken into account following the lines of by introducing the real 
part of the pion selfenergy in the pion propagators cut by the horizontal line 
in Fig. ^ when one evaluates the photon selfenergy with a ph and two pions 
in the intermediate state. We find a moderate increase of the cross section 
by an amount of 20 percent in ^^C from this effect when the tttt interaction 
in the medium is considered in addition, but again it does not modify at all 
the shape of the invariant mass distribution. All these things discussed, we 
can conclude that while there are uncertainties of about 20 percent in ^^C, 
and a bit more in ^°^P6, in the size of the total cross section, the shape of the 
mass distribution is still determined basically by the implementation of the 
TTTT interaction in the medium. 

In Fig. 1^ we can see the results for the invariant mass distribution of the 
two pions for and ^osp^ including the absorption of the final pions. The 
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Figure 7: Two pion invariant mass distribution for 2tt^ photoproduction in ^^C and 
"^^^Pb. Continuous lines: Using the in medium final vrvr interaction. Dashed lines: 
using the final tttt interaction in free space. 



difference between the solid and dashed curves is the use of the in medium vrvr 



scattering and G function instead of the free ones, which we take from |17|. 
As one can see in the figure, there is an appreciable shift of strength to the 
low invariant mass region due to the in medium vrvr interaction. This shift is 



remarkably similar to the one found in the preliminary measurements of [ff? 

These results show a clear signature of the modified vrvr interaction in the 
medium. The fact that the photons are not distorted has certainly an advan- 
tage over the pion induced reactions and allows one to see inner parts of the 
nucleus. In this sense it is worth noting that from our calculations we can 
determine the average nuclear density felt by the reaction which turns out to 
be 35 per cent and 65 per cent of the normal nuclear density for ^^C and "^^^Pb 
respectively. 

Although we have been discussing the nir interaction in the nuclear medium 
it is clear that we can relate it to the modification of the a in the medium. 
We have mentioned that the reason for the shift of strength to lower invariant 
masses in the mass distribution is due to the accumulated strength in the 
scalar isoscalar vrvr amplitude in the medium. Yet, this strength is mostly 
governed by the presence of the a pole and there have been works suggesting 
that the a should move to smaller masses and widths when embedded in the 



nucleus gl], |8l g|, |3g. We should stress that in refs. |31|, g|, |2 
effects of chiral symmetry restoration from the dropping of the condensate 
and/or the change in the bare a mass are considered. In our approach which, 
as we mentioned, relies upon a conventional many body expansion based on 
standard chiral lagrangians, where the a in the free space or in the nucleus is 
generated dynamically, we can look at the a poles and see their evolution with 
the nuclear density. A detailed study of the poles in the complex plane for the 
TTTT interaction in the nuclear medium within this approach is now available 
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in 



50|| . In this work it is indeed found that the pole position of the a moves 
to smaller energies as the density increases and the width is also reduced. 
The present results, when contrasted by the definitive data, if they confirm 



the preliminary ones of |T^, should represent an evidence of this interesting 



phenomenon which would come to strengthen once more the nature of the 
a meson as dynamically generated by the multiple scattering of the pions 
through the underlying chiral dynamics. 
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